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Abstract 

In this paper we consider a reaction-diffusion equation of Fisher-KPP type 
inside an infinite cylindrical domain in coupled with a reaction-diffusion 

equation on the boundary of the domain, where potentially fast diffusion is 
allowed. We will study the existence of an asymptotic speed of propagation for 
solutions of the Cauchy problem associated with such system, as well as the 
dependence of this speed on the diffusivity at the boundary and the amplitude 
of the cylinder. 

When N = 1 the domain reduces to a strip between two straight lines. 

This models the effect of two roads with fast diffusion on a strip-shaped field 
bounded by them. 
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1 Introduction 


Recently, in |3l IH |5], the authors introduced a model to describe the effect of a road 
with potentially fast diffusion in a field which was assumed to be the (half) plane. 
The scope of the present paper is to provide a rigorous mathematical framework for 
the problem of the propagation of biological species or substances in presence of more 
general domains, with potentially fast diffusion on their boundary. 

In particular we will consider the case of diffusion inside a {N + 1)-dimensional 
cylinder = M xBn{0, R), where Bn^O, R) is the A^-dimensional ball of radius R > 0 
centered at the origin. In this case the walls, which may favor the process of diffusion, 
play the role of a pipe. 

In the specific case = 1, this situation models two parallel roads of potentially 
fast diffusion bounding a field described be a two-dimensional slab = M x (—R, R). 
Such spreading heterogeneities have been tracked in several real environments, such 
as the diffusion of wolves along seismic lines in the Western Canadian Forest m) 
or the early spreading of HIV among humans in Democratic Republic of Congo ([3)) 
just to mention a few ones (see [3] for a more detailed list of empirical observations). 

In our model, roughly speaking, a species, whose density will be denoted by v, 
occupies the region D, where it reproduces, diffuses and dies, according to a logistic 
law. The logistic law is governed by a function / (as a model case, the reader may 
think about the case of f{v) = u(l — u), but more general logistic laws will be allowed 
as well in this paper). 

The population diffuses randomly inside the field, with a diffusion coefficient d. 
When the individuals of the species, driven by their random movement, hit the bound¬ 
ary of D, a proportion u of them decides to use it to keep diffusing, with a possibly 
different diffusion coefficient D. We denote the density of the population which moves 
on dQ by u and we assume that there no reproduction occurs, but the number of 
individuals on it only varies because of the above mentioned positive contribution 
from the field and a negative contribution given by a part /i of the population that 
leaves dQ to go back to the field. 

Motivated by the wide number of applications, several authors have recently con¬ 
sidered, from different point of views, similar models which couple reaction-diffusion 
equations posed in different spatial dimensions and with different diffusivities in the 
interior and on the boundary of a bounded domain. For example, in |2], the authors 
consider no reaction in the interior and nonlinear Robin coupling terms. They prove 
the existence of a unique weak solution of the evolution problem by means of sub- 
and supersolutions and determine the exponential rate of convergence to the unique 
steady state through entropy methods. On the other hand, in [10], the authors con¬ 
sider two coupled reaction-diffusion systems, one in the interior of the domain and 
one on the boundary and study, according to the values of the different diffusivities, 
the formation of Turing patterns in the interior and/or the boundary. 

Our first result shows, roughly speaking, that the species invades the whole en¬ 
vironment O, i.e. it converges to a positive steady state of our model (which will be 
shown to be unique). With this result in hand, a very natural question from the bio¬ 
logical point of view is to know how fast the species invades the domain, or, in other 
words, what is the velocity of expansion of the region in which the solution of our 
problem converges to the positive steady state. Such concept is known as asymptotic 
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speed of propagation and it was introduced in the contest of population dynamics 

in pin]. 

The main results of this paper are that our problem indeed admits a positive 
asymptotic speed of propagation c* and we study its dependence on the parameters 
d, D and R, characterizing the different regimes of propagation. In particular, we will 
prove that, taking a large diffusion coefficient D produces arbitrarily large speeds, but 
taking D small cannot slow down the process below a certain threshold. Notice that 
neither of these two phenomena is obvious a priori from the biological interpretation 
of the model, since, in principle, one can think that the slower diffusive process might 
drag the whole diffusion and either prevent arbitrarily large speeds of propagation or 
make such speed arbitrarily small - but we will show that this is not what happens 
in our model. 

Moreover, we will prove that increasing the diffusion coefficient D always makes c* 
larger, as illustrated in Figure [T] Remarkably, this phenomenon occurs even if D is 
very small, in particular smaller than the diffusion coefficient d inside the held (i.e. the 
speed of propagation is inhuenced by the auxiliary boundary device even when the 
latter is not, in principle, faster than the held itself). This is in contrast with the 
case of a single road treated in |3|, where, for D < 2d, the speed of propagation is 
constantly equal to ckpp, the asymptotic speed of a classical Fisher-KPP equation in 



Figure 1: Behavior of the asymptotic speed of propagation c* with respect to D, the 
diffusivity on the boundary. 


The dependence of c* on the geometry of the held is even more intriguing and rich 
of somehow unexpected features. If the held is very thin (i.e. when R goes to zero), 
the asymptotic speed approaches zero. Such fact is not immediately evident from 
the biological interpretation: in the case of two fast dihusive roads, for instance, one 
might be inclined to think that for small R the two roads get very close and more 
“available” to the individuals, which should make the dihusion faster. This is not the 
case: the explanation of this phenomenon probably lies in the fact that small helds 
do not allow the species to reproduce enough, and this reduces the overall diffusion. 
On the other hand, as R becomes large, the asymptotic speed c* approaches the limit 
value c^, the asymptotic speed in the half-plane found in [3|. As mentioned above, 
coincides with the standard speed ckpp in the held if D < 2d. Otherwise, if D > 2d 
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(i.e. when the auxiliary network provides fast enough diffusion, with respect to the 
diffusion in the held) is greater than ckpp- 

The threshold D = 2dis also crucial for the monotonicity properties of the asymp¬ 
totic speed c* in dependence of the size of the held R. Indeed, when D < 2d we have 
that c* is increasing in R (see the left graph in Figure]^. This can be interpreted 
by saying that, if the dihusion on the auxiliary network is not fast enough, then the 
propagation speed only relies on the proliferation of the species, which is in turn 
facilitated by a larger held. Conversely, when D > 2d, the speed c* does not depend 
monotonically on R anymore, but it has a single maximum at a precise value R = Rm, 
as shown in the right graph of Figure 




Figure 2: Behavior of c* with respect to the amplitude of the cylinder in the case D < 2d 
(left) and D > 2d (right). 


Once again, this fact underlines an interesting biological phenomenon: when the 
auxiliary network provides a fast enough dihusion, there is a competition between the 
proliferation of the species and the availability of the network, in order to speed up the 
invasion of the species. On the one hand, a large R would favor proliferation and so, 
in principle, dihusion; on the other hand, when R becomes too large the fast dihusion 
device gets too far and cannot be conveniently reached by the population (similarly, 
for small R the fast dihusion device becomes easily available, but the proliferation 
rate gets reduced). As a consequence of these opposite ehects, there is an optimal 
distance at which the two roads shall be placed to get the maximum enhancement of 
the propagation. 

Such kind of monotonicity might arise also in other contests with contrasting 
ehects on the propagation, for example if D > d and there is a smaller reaction 
on the boundary of the cylinder than in the interior, or if D < d and the reaction 
on the boundary is larger. However, this kind of analysis, as well as the ehect of 
transport terms on the boundary, goes out the scope of this work and is left for 
further investigation. 

After these considerations, we pass now to introduce the formal mathematical 
setting that describes our model and state in detail our main results. 

We consider the following problem 

(vt{x,y,t) - dAv(x,y,t) = f(v(x,y,t)) for (x, ?/) G H, t>0 

<Ut{x,y,t)-DAoQu{x,y,t) = iyv{x,y,t)-pu{x,y,t) for {x,y) e dn, t>0 (1.1) 

{ddnv(x,y,t) = pu(x,y,t) — iyv(x,y,t) for (x,y) G dfl, t > 0, 

where denotes the Laplace-Beltrami operator on dQ, while dn stands for the 
outer normal derivative to dfl. The parameters D,d, yi,v,R are positive and the 
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nonlinearity / G C^([0,1]) is of KPP type, i.e. satisfies 

/(0) = 0 = /(l), 0 </(s) </'(0)s for se (0,1) (1.2) 

and is extended to a negative function outside [0,1]. When iV = 1, due to the non¬ 


connectedness of 9i?i(0,i?), problem (1.1) reduces to the case of a strip bounded by 
two roads and reads as 


'vt{x, y, t) - dAv{x, y, t) = f{v{x, y, t)) for (x, y) G1 
Ut{x, t) — Duxx{x, t) = uv{x, R, t) —y,u{x, t) for a; G M, 

dvy{x, R, t) = yu^x, t) — uvi^x, R, t) for a; G M, 

Ut{x, t) — Duxx{x, t) = i'v{x, —R, t) —yu{x, t) for a; G M, 

—dvy{x, —R, t) = yui^x, t) — uv{x, —R, t) for a: G M, 


.x{—R, R), t>0 
t>0 
t>0 
t>0 
t>0. 


(1.3) 


This kind of models goes back to [3], where the authors studied the problem 

^vt{x,y,t) - dA^N+iv{x,y,t) = f{v{x,y,t)) for (a:, ?/) gM^ x IR+, f > 0 
Ut{x, t)—DA]^Nu{x, t) = iyv{x, 0, t) —yu{x, t) for x G f > 0 

[ddnv{x, 0, t) = yu{x, t) — uv{x, 0, t) for x G 


iN 


(1.4) 


f > 0 


for = 1, which is the counterpart of (1.3) with just one road and the field which 


is a half-plane. Some results about (1.4 


, such as the well-posedness of the Cauchy 
douville-type result for stationary solutions. 


problem, the comparison principle and a 
have been proved in |3] in arbitrary dimension N. 

The main result of |3] (proved there for = 1, but the arguments extend to 
arbitrary dimension, as shown in Section below) is that system (1.4) admits an 
asymptotic speed of spreading c^, along any parallel direction to the hyperplane 
{?/ = 0}. If H < 2d, coincides with ckpp, where 


ckpp 


= 2 ^/^) 


is the spreading speed for the Fisher-KPP equation Vt—dAv = f{v) in the whole space 
]^Af+i j-ggg while, if H > 2d, is strictly greater than cxpp. As a consequence, 
the results of |3] show that a large diffusivity on the road enhances the speed of 
propagation of the species in the whole domain. 

The main results of the present paper concern the asymptotic speed of spreading 


for problem (1.1). The first one is related to its existence. 


Theorem 1.1. There exists c* = c*{D,d, y,^, R, N) > 0 such that, for any solution 


{u,v) of (1.1) with a continuous, nonnegative, initial datum {uo,vo) ^ (0,0), with 
bounded support, there holds 


(i) for all c> c* 


lim sup \u{x,y,t)\ = 0, 

\x\>ct 

\y\=R 


lim sup \v{x,y,t)\=0, (1.5) 

\x\>ct 

\y\<R 


(a) while, for all 0 < c < c*, 


lim sup 

|x|<ct 

\y\=R 


u{x,y,t) — 
y 


= 0, lim sup \v{x,y,t) — 1\ = 0. (1.6) 

|x|<ct 

hl<« 
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Next, we derive the following qualitative properties for the speed of spreading. 


Theorem 1.2. The quantity c* given by Theorem (1.1) satisfies the following prop¬ 
erties: 

(i) for fixed d, /r, z/, R, N, the function D i—>• c*{D) = c*{D, d, /r, z/, R, N) is increas¬ 
ing and satisfies 


lim c*(Zl) =: Cq > 0, 

DiO ^ 


, c*(D) , 

lim —G (0 ,+cxd). 


(1.7) 


D—^+oo 

In particular c*{D) —>■ +cx) as D ^ +cx); 

(a) for fixed D, d, /i, u, N, the function R i—)■ c*{R) = c*{D, d, fi, z/, R, N) satisfies 

( 1 . 8 ) 


hnic*(i?) = 0, 

Rl-O 


lim c*{R) = c^, 

R-^+oo ^ ^ °° 


where is the asymptotic speed of spreading for problem (1.4), which satisfies 

= ckpp if D < 2d, 

> Ckpp if D > 2d. 

Moreover, if D < 2d, the function R ^ c* (R) is always increasing, while, if 
D > 2d, it is increasing up to the value 

NDu 


Rm 


{D — 2(i)/i’ 

and decreasing for greater values. In addition, 

c*(Rm) = maxc*(R) = . = 

^ ^ «>o ^ ^ ^Ad{D - d) 


Ckpp- 


(1.9) 


( 1 . 10 ) 


The paper is organized as follows: in Section we present some preliminary 
results related to the existence and comparison properties of solutions of the Cauchy 
problem associated with (1.1). In Section]^ we prove the existence and uniqueness 
of the positive steady state and we show that it is a global attractor, in the sense 
of locally uniform convergence, for solutions of (1.1) starting from a nonnegative, 
nontrivial initial datum. Section is dedicated to the geometrical construction of c* 
and to the proof of Theorem |1.1[ Finally, the dependences of c* on D and on R are 
discussed in Sections and respectively. 


2 Preliminary results and comparison principles 

In this section we present some fundamental results that are extensively used in the 
rest of the paper. Some of them are contained in or follow easily from [31 H in 
which cases we refer to such works and just outline the adaptations that need to be 
done. We start with a weak and strong comparison principle for the following system 

{ Vt — dAv + cvx = f{v) in 12 X (0, +cx)) 

Ut — DAqqu + cUx = uv — p,u in dH x (0,+c5o) (2.1) 

ddnV = p,u — uv on dTl x (0, +cx)). 
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which is more general than (1.1), since we include the possibility of an additional 
drift term with velocity c G M in the x-direction. As usual, by a supersolution (resp. 


subsolution) of (2.1) we mean a pair {u,v) satisfying System (2.1) with “>” (resp. 


“<”) instead of “=”. Below and in the sequel, the order between real vectors is 
understood componentwise. 

Proposition 2.1. Let {u,v) and {u,v) be, respectively, a subsolution bounded from 


above and a supersolution bounded from below of (2.1) satisfying {u,v) < {u,v) at 
t = 0. Then {u,y) < {u,v) for all t>0. 

Moreover, if there exists T > 0 and {x, y) E fl such that either u{x, T) = u{x, T) 
or yfx, y, T) = v{x, y, T), then {u, v) = (u, v) for t G [0, T]. 

Proof. The arguments of the proof of |3], Proposition 3.2] easily extend to our case, the 
key point being the strong maximum principle, which holds for the Laplace-Beltrami 
operator Ago (and the associated evolution operator) like for the operator d^x- □ 

The following result establishes the existence and uniqueness for the solution of 


the Cauchy problem associated with (1.1). The proof of the existence part follows 


from [31 Appendix A], while the uniqueness follows from Proposition 2.1 


Proposition 2.2. Let {uq{x, y), Vq{x, y)) be a nonnegative, bounded and continuous 
pair of functions defined in dfl and 12 respectively. Then, there is a unigue nonnega¬ 


tive, bounded solution {u,v) of (2.1) satisfying {u,v)\t=o = {uo,vo). 


Once these properties established, we can show that the model exhibits total mass 
conservation if no reproduction is present in the field. 


Proposition 2.3. Assume / = 0 and let {u,v) be the solution of (1.1) with a non¬ 
negative, bounded initial datum (uo,Vo) decaying at least exponentially as |a;| —)■ oo. 
Then, for every t > 0, we have 




Li(0) 


+ 


Li(ao) “ Ikoll Li(0) + ll«o|| LI (90) ' 


Proof. On the one hand, taking as subsolution the pair (0,0) in Proposition 2.1 we 
see that {u, v) is nonnegative for all times. On the other hand, anticipating on Section 
the system admit exponential supersolutions with arbitrary slow exponential decay 
(i.e., with the notation of Section]^ a arbitrarily small, which is allowed provided c is 
sufficiently large) and therefore, up to translation, above (mq, Vq) at t = 0. Proposition 
then implies that {u,v) decays exponentially as |a;| —)■ oo for all times. As a 


2.1 


consequence of standard parabolic estimates, the same is true for the derivatives of 
u and V, up to the second order in space and first order in time. Therefore, using the 


equations of (1.1) and applying Stokes’ theorem, we obtain 


d 

dt 


^ll^^(^^)llLl(o)= / vt{x,y,t) = d Av{x,y,t) = d dnv{x,y,t) 

^ ^ JQ JQ Jan 

(nu{x,y,t)-uv{x,y,t)) = D A3nu{x,y,t)- Ut{x,y,t) 
f Jan Jan 


Jan 

^)llLi(aQ) 
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where, letting Aqn-i denote the Laplace-Beltrami operator on the {N — l)-dimen- 
sional unit sphere, we have used 


r r r-\-oo p-\-oo r 

/ AQQu{x,y,t)= / dl^u{x,y,t)+ / AsN-iu{x,y,t) = 0, 
Jdn J-OO 


which holds because d^u decays exponentially to 0 and the second addend in the right 
hand side is 0 by Stokes’ theorem, since is compact and has no boundary. We 

have eventually shown that is constant in time and this 

concludes the proof. □ 


We also need the following comparison principle involving an extended class of 
generalized subsolutions and which is a particular instance of [5|, Proposition 2.2]. 
Actually the result of [5] holds for a more general class of subsolutions than the one 
we consider here, however we present it in the form needed in the sequel. 


Proposition 2.4. Let (Mi,ni) he a subsolution of (2.1), bounded from above and such 
that Ui and Vi vanish on the boundary respectively of an open set E of dfl x [0, +oo) 
and of an open set F of fl x [0, +cxd) (in the relative topologies). If the functions u, 
V defined by 


u 


max{Mi, 0} 
0 


in E 

otherwise, 


V 


maxjni, 0} in F 
0 otherwise, 


satisfy 


v{x,y,t) > vi{x,y,t) 
uix,y,t) > ui{x,y,t) 


for all {x,y,t) G dQ x M_|_ such that u{x,y,t) > 0, 
for all {x,y,t) G dfl x IR+ such that y{x,y,t) > 0, 


( 2 . 2 ) 


then, for any supersolution {u,v) of (1.1) bounded from below and such that {u,y) < 
{u, v) att = 0, we have {u, y) < (u, v) for all t > 0. 


Before concluding this section, we derive the strong maximum principle for two 
variants of system (1.1): 


(vt- 

- dAv = 

= fix) 

for 

ix, 

y) e n, 

t 

> 0 


\ut -- 

= uv — _ 

yu 

for 

ix, 

y) e dfl, 

t 

> 0 

(2.3) 

[dd„ 

,V = piU 

— uv 

for 

ix, 

y) e dfl, 

t 

>0, 


(vt- 

- dAyV 

= fix) 


for 

{x,y)efl, 


f>0 


ut- 

-AqqU -- 

= z/n — , 

yu 

for 

ix,y)e dfl, 

o 

A 

(2.4) 

[ddn 

V = pLU 

— i/v 


for 

(x,2/)g dfl, 

o' 

A 



where Ay denotes the Laplace operator with respect to the ^/-variables only. Both 
systems are semi-degenerate, in the sense that exactly one of the two parabolic equa¬ 
tions is degenerate: the second one for (2.3) and the first one for (2.4). Apart from 
being in themselves interesting, these results will be needed in Section to charac¬ 
terize some singular limits related to the asymptotic speed of propagation for (1.1). 
Roughly speaking, even if one equation is degenerate, the presence of the other, which 
is not degenerate, still guarantees the system to be strongly monotone. 








Proposition 2.5. Let {u,v) and {u,v) be a subsolution bounded from above and a 
supersolution bounded from below, both of system (2.3) or of system (2.4), such that 
{u,v) < (u,v) in fl X (0, +cxd). If there exist T > 0 and {x,y) G 12 for which either 
u{x, T) = u{x, T) or yfx, y, T) = v{x, y, T), then {u, v) = (u, v) for t G [0, T]. 

Proof. The pair {u, v) := {u — u,v — v) is a nonnegative supersolution of either system 


(2.3) or (2.4), with f{v) replaced by v{f{v) — f{y))/{v — y) in the hrst equation. This 


replacement will be understood throughout the proof. We argue slightly differently 
according to which system is involved. 

Case 1. System (2.3). 

Suppose hrst that v{x,y,T) = 0. If {x,y) G 12 then, applying the strong parabolic 


maximum principle to the hrst equation in (2.3), we conclude that f = 0 for t <T 
and thus the same is true for u by the third equation. If {x, y) G c212, the parabolic 
Hopf lemma implies that either dnv{x,y,T) < 0 or n = 0 for t < T. The hrst 


situation being ruled out by the third equation in (2.3), we can conclude as before. 


Suppose now that u{x,y,T) = 0. We claim that this yields v{x,y,T) = 0, and thus 
the previous arguments apply. Indeed, if it were not the case, we would have from 


the second equation of (2.3) that Ut{x,y,T) > 0 and, as a consequence, u{x,y,t) < 0 
for 2 < T, 2 ~ 0, which is impossible. 


Case 2. System (2.4). 


Suppose hrst that u{x,y,T) = 0. By the second equation in (2.4) we see that u 


is a supersolution of the uniformly parabolic linear operator dt — Aqq + p. Hence, 
the parabolic strong maximum principle yields m = 0 in 2212 x [0,T], and thus the 


same is true for v, always by the second equation in (2.4). Moreover, dnV = 0 on 
2212 X [0,T] by the third equation in (2.4). Since, for hxed a: G M, {y,t) h- >■ v{x,y,t) 
is a supersolution of the hrst equation in (2.4), which is uniformly parabolic in {y,t), 
we deduce from the Hopf lemma that n = 0 in 12 x |0,T]. 

Consider now the case where v{x,y,T) = 0. If {x,y) G 2212, since n > 0, we have 


that dnv{x,y,T) < 0 and the third equation of (2.4) implies u{x,y,T) = 0. That is. 


we end up with the previous case. If instead (£, G 12 then, applying the parabolic 
strong maximum principle to {y,t) i—>■ v{x,y,t), which is a supersolution of the hrst 


equation in (2.4), we derive v{x,y,t) = 0 for \y\ < R, t < T, i.e. v vanishes on a 
point of 2212, which is a case that we have already treated. □ 


Remark 2.6. By looking at the proof of Proposition |2.5[ it follows immediately that 
the comparison principles also hold if we add to ( |2.3| ) and (2.4) a constant drift term 
in the x-direction. 


3 Liouville-type result and long time behavior 


In this section we discuss the asymptotic behavior as 2 —)■ +oo of solutions of the 
Cauchy problem. First of all, for later purposes, we study the long time behavior 


of solutions of Problem (1.1) with an additional drift term and which start from a 
specihc class of initial data. Next, we show that the limit in time of a solution of 


(1.1) is always constrained between two positive steady states that do not depend on 


the x-variable and are rotationally invariant in y. Finally, we classify all the steady 
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states with these symmetries, showing that there is a unique nontrivial one, which is 
a global attractor for our problem. 


Proposition 3.1. Let {u,v) ^ (0,0) be a nonnegative generalized stationary subso¬ 


lution, in the sense of Proposition 2-4, of (2.1) which has bounded support and is 


rotationally invariant in y. Then, the solution {u,v) of (2.1) with {u,v) as initial 
datum satisfies 

hm {u,v) = {U,V), (3.1) 

C^ + OO 


locally uniformly in 12, where {U,V) is a positive stationary solution of (2.1), which 
is rotationally invariant in y and x-independent. 


Proof. By the comparison principle of Proposition 2.4 
Hence, by Proposition 2.1, {u,v 


(m, v) > (m, v) for all times. 


is nondecreasing in time, and it is actually strictly 
increasing by the strong comparison result, because otherwise {u,v) = {u,y) for all 
times, which is impossible since (m,I t) is not a solution of (1.1). From this mono¬ 


tonicity in t, we obtain that {u, v) tends, as t —>■ -foo, locally uniformly in 12 to a 


nonnegative bounded steady state (U,V) of (2.1), which proves (3.1). 

As for the symmetries of (27,17), we have that it is rotationally invariant in y since 
the Cauchy problem itself is rotationally invariant in y and has a unique solution. To 
prove that it is also ^-independent, we “slide” the initial datum {u,y) as follows. 

From the monotone convergence, we have that (27,17) > (m, v) and, since the latter 
has compact support, (M(-±h, ■),i;(-±h, ■)) still lies below (27,17) for sufficiently small 
h. It follows from Propositions 2.4 and 2.1| that the solution to (1.1) emerging from 
± ± h, ■)) lies below (27,17) for all times. But, by the x-translation 

invariance of the system, such solution is simply the translation in the x variable by 
±h of (m, v) and therefore it converges to the corresponding translation of (27,17) as 
t —?• -|-oo. We have eventually shown that these translations of (27,17) he, for h small 
enough, below (27,17) itself, which immediately implies that (27,17) does not depend 


on X. 


□ 


Theorem 3.2. For every nonnegative, bounded (uo,Vo) ^ (0,0), there exist two pos¬ 


itive steady states of (1.1), denoted by (27i, fo) and (272, lA); which are x-independent 


and rotationally invariant in y and such that the solution of (1.1) starting from the 
initial datum (Mo,fo) satisfies 


(27i,17i) < liminf(M, n) < hmsup(M,n) < (272,172), 

t^-^OO t^ + oo 


(3.2) 


where the first ineguality holds locally uniformly in VL and the last one uniformly. 
Proof. It is easy so see that the pair {u,v), with 


V = max <jl, iluoiloo, ^il^oilc 


_ z/_ 

u = —n, 

/i 


is a stationary supersolution to (1.1) which is above {u, n) at t = 0. As a consequence 


of Proposition 2.1, we have that the solution of (1.1) starting from {u,v) lies always 


above {u,v) and tends non-increasingly, as t —)■ -|-cx), to a positive bounded steady 


state (272, 1 ^ 2 ) of (1.1). It also follows from the symmetries of the system, that (272,14) 
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inherits from (m, v) the x-independence and the radial symmetry in p. As a conse¬ 


quence, the convergence is uniform in Q and the last inequality in (3.2) is thereby 
proved. 


To prove the hrst one, consider such that the eigenvalue problem 


— A(j) = (j) in i?7v(0, i?i), 

0 = 0 on 9i?Ar(0, i?i). 


(3.3) 


admits a positive solution 0. Such eigenfunction is therefore the (unique up to a 
scalar multiple) principal eigenfunction of —A in i?Ar(0, Ri) under Dirichlet boundary 
condition, and the associated eigenvalue is 1. It is well known that 0(j/) = 'ip{\y\) with 
0 : (0,i?i) —)■ M is decreasing and satishes 0'(O) = 0. Consider now the pair {u,v) 
dehned by 

/ cos(ax) ( 1 , 70 ( 01 /)) if |x| < ^ 

~ ' \ ( 0 , 0 ) otherwise, 

with a, 7 > 0 and 0 < 0 < Ri/R to be chosen. Plugging (m, v) into the third equation 
of (1.1), we obtain 

7 = - - - (3.4) 

^ d 00 '( 0 i?)z/ 0 ( 0 i?) ^ ^ 

which is positive for 0 < 0 < 0 < where 0 is the first positive zero of the 

denominator. We now look for a, 0 so that {u, v) is a generalized subsolution - in the 


sense of Proposition 2.4 - to 


Vi — dAv = 

Ut - DAqqU = - l^u. 


Due to (3.4), this results in the system 

da^ + d/3^ < /'( 0)/2 


DW < 


—pd(3'ilj\l3R) 

dl3'ifj'{/3R) + i"ifj{/3R) 


(3.5) 


For fixed 0 < min 


/'(o) 

2d • 


, since the right-hand side of the second inequality of 


(3.5) is positive, it is possible to take a ~ 0 in such a way (3.5) is satisfied. Thanks to 


( 1 . 2 ) we have that, for e ~ 0 , e{u, v) is a compactly supported generalized subsolution 
to (1.3). Reducing e if need be, we can further assume that e{u,v) lies below the 
pair {u,v) at time 1, the latter being strictly positive by Proposition 2.1 Again 


by Proposition 2.1, the order is preserved between {u,v) shifted by 1 in time and 
the solution {u,v) to (1.1) emerging from e{u,v). The proof is concluded applying 
Proposition 3.1[ 


□ 


The previous result indicates that we have to focus the attention on the steady 


states of ( 1 . 1 ) or, more generally of ( 2 . 1 ), with the symmetry properties specihed 


there. In this sense we have the following Liouville-type result. 


Proposition 3.3. The unique nonnegative bounded stationary states of (2.1) which 
are x-independent and rotationally invariant in y, are ( 0 , 0 ) and ( 751 )- 
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Proof. Let ( m , v) be a steady state as in the statement of the proposition. Then u 
is constant and there exists a nonnegative function T G C^([0,/2]) nC^((0,i?]), with 
'L'(O) = 0, such that v{x,y) = 'L(||/|). It follows from the second equation in (2.1) 


that u = and thus the other two equations read 


-d^"^r) - = /(T(r)) 

T'(i?) = 0 = T'(0). 


r e (0,i?), 


(3.6) 


To prove the result it is sufficient to show that = 1 or T(i?) = 0. Suppose that 

this is not the case. Then ^ 0 and T' does not vanish in a left neighborhood 

of i?, being positive if 0 < < 1 and negative if > 1. Set 

p := max{r G [0, R) : \I/'(r) = 0}. 

From this dehnition we have that T' has a hxed strict sign in (p, R), which is the same 
as 1 — If this sign is positive then, for p < r < R, we have that 0 < T(r) < 

T(i?) < 1, whence /(^(r)) > 0 and the hrst equation in (3.6) eventually yields 


T"(r) < 0. This is impossible because t&'(p) = = 0. If instead 1 — '^{R) < 0, 

then we obtain T'(r) < 0 for p < r < /?, which implies T(r) > $(«) > 1 and thus 


/(T(r)) < 0. Therefore, in such case, the hrst equation in (3.6) yields T"(r) > 0 for 
p < r < R, which again contradicts tb^p) = T'(i?) = 0. □ 


As an immediate consequence of Theorem 3.2 and Proposition [3^ we characterize 
the long time behavior of solutions of the Cauchy problem associated with (1.1). 


Corollary 3.4. Any solution {u,v) of (1.1) starting from a bounded, nonnegative 
initial datum {uq,vo) ^ (0,0) satisfies 


V 


lim u{x,yR) = -, 

t —^ + oo 

lim v(x, y, t) = 1, 

t^+oo ^ ^ ’ 


locally uniformly in dVl, 
locally uniformly in VL. 


(3.7) 


In particular, restricting to stationary solutions of (1.1), we obtain that the unigue 


nonnegative bounded steady states are (0,0) and 1 


Remark 3.5. (i) Observe that this last statement of Corollary |3.4| is much stronger 

than Proposition |3.3 with c = 0, because it holds without knowing a priori the 
symmetry of solutions. 


ii) In (3.7), thanks to Theorem 3.2, the inequalities “<”related to the limsup hold 


uniformly in O. 


4 Asymptotic speed of spreading 


In this section we prove the existence of a value c* such that (1.1) admits a super¬ 


solution moving with speed c* and some generalized subsolutions moving with speed 
less than and arbitrarily close to c*. This value c* will be identihed as the asymp¬ 
totic speed of spreading appearing in Theorem 1.1 The construction of c* will also 
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provide some key information about its dependence on D and R that will be used in 


the following sections to derive Theorem |1.2 

The starting point to hnd c* is the analysis of plane wave solutions for the lin¬ 
earization of (1.1) around n = 0. This is achieved in Section [4T| through a geometrical 
construction. The plane waves are supersolutions of (|1.1[ ) by the KPP hypothesis and 
their existence immediately implies (1.5). In Section 4.2 we construct the generalized 


subsolutions and prove ( 1 . 6 ) 


4.1 Plane wave solutions 


Consider the linearization of (1.1) around n = 0: 


Vt — dAv = f'{0)v in X M 

Ut — DAgfiu = vv — fiu on dVl x 
ddnV = jjiu — vv on dVt x 

We look for plane wave solutions in the form 

{%v) = e"(^+^*^(l, 70 (/?,|/)), 


(4.1) 


(4.2) 


with a, 7 > 0 , /3 G M and 0 positive, i.e., moving leftward at a velocity c and decaying 
exponentially a.s x + ct —?• —oo. In contrast with [3], but in analogy with [13], we need 
to consider two types of plane waves, corresponding to a dichotomy in the dehnition 
of 0 : 

[01 (01/) if0>O 


0 ( 0 , 2 /) : = 


[ 02 ( 02 /) if 0 < 0 . 


The functions 0i and 02 are related to the eigenvalue problem for the Laplace operator 
in a ball and in the whole space respectively. We will show that 0 is differentiable, 
reflecting a continuous transition from 0 i( 0 //) to 02(02/) at 0 = 0 , uniformly in y, 
due to the fact that the support of 0i(02/) becomes the whole as 0 0 0. The 
construction of the subsolutions that will be carried out in Section 14.21 makes use of 
one or the other type of plane wave depending on the values of the parameters of 
the problem. This dichotomy will give rise to the two different monotonicities with 


respect to R stated in Theorem 1.2 ii). 

We start with 0 > 0. In this case we take (j){(3,y) := 0i(02/) in (4.2), where 


01 is the positive eigenfunction 0i of problem (3.3), normalized by 0i(O) = 1. We 
therefore impose 0 < Ri/R. We recall that (j)i{z) = 'ipi{\z\) for a real analytic 
decreasing function 0i on [0,i?i] such that 0'i(O) = 0. Plugging this expression into 


(4.1) we are driven to the system 


—da'^ + (i0^ - 1 - ca = /'(O), 
—Da^ - 1 - ca = 1^701 (0i?) — /x, 
d'^l3il)'^{ldR) = /i — 1^701 (0i?). 


(4.3) 


Solving the last equation for 7 yields 

7 = 




(i00'l(0i?) -b 1X01 (0i?) 


(4.4) 
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In order to have 7 > 0 we restrict to 0 < /3 < /?, where f3 G (0,-Ri/-R) is the first 
positive zero of the function r 1 —>■ dr'ijj[{rR) + which is positive for r = 0 

and negative for r = Ri/R. By plugging (4.4) into the second equation of (4.3) we 
obtain 

—fid(3'ilj[{j3R) 


-Da + ca = 


d(3'ip[{(3R) + h''ipi{(3R) 


=■ XiW) 


and, solving for a, 


a 


(c, ^) = ^ (‘^ =*= 7 ^ - 4C\'i 


(4.5) 


Observe that Xi is positive in (0,/9), satishes Xi(0) = 0 = x'i(O) and XiW) +cx) as 
/3/3. We now use a property of the function that will be crucial also in the sequel: 
log'll is concave (see e.g. 0 ). It implies that the negative function 'ipi{r)/'ip[{r) is 
increasing and then the same is true for r 1 —>■ 'ipi{r)/. Reasoning on 1/xi, we 
hnd that Xi is increasing too. Therefore, for every c > 0 there exists a unique value 
of G (0,/d), denoted by /3{c), such that 


= 4:Dxi0{c)). 


There holds 


lim/3(c) = 0 , 
c4,o 


lim /3(c) = /d. 

C-/. + 00 ^ ^ 


The functions are real-valued for /d G (0,/3(c)], where they satisfy 

(3 H->• /3) is decreasing, [3 1 —>■ a]/)(^c,/3) is increasing. 

With regard to the monotonicity in c, we have that 

c H->■ a'^/c,/3) is increasing, c h-)■ a]/){c,/3) is decreasing 

and that the region delimited by q:^(c, •) invades, increasingly, the strip [0,/S] x 
(0, -|-cx)) in the (/S, Q:)-plane as c —>■ -|-cx), while it shrinks to (0,0) as c 0. 

On the other hand the hrst equation in (4.3) represents the two branches of a 
hyperbola 


(c, /3) := ^ (c ± - 4 pp + 4 ^ 2/52 j 


where ckpp = 2^Jdf'{0). For c < ckpp, the functions a^{c,/3) are real-valued for 
/S > /S(c), where /S(c) >0 satishes 


c = c 


KPP 


4(i^/3(c)^. 


Observe that /S(c) —>■ f'{0)/d as c 0. On the contrary, for c > ckpp, the hyperbolas 
are dehned for every /3 > 0. In particular, for c = ckpp, the hyperbolas degenerate 
into the straight lines ±/3 -t-CKpp/( 2 d). As c increases to -|-cx), the region lying in the 
hrst quadrant between the curves a'^(c, •) with /3 > /3(c) ii c < ckpp and /S > 0 if 
c > Ckpp, invades monotonically the whole quadrant. 

Passing to the construction of the second type of plane waves, we consider the 


positive radial eigenfunction 02 of A 02 = 


m 


t,N 


, normalized by 02 ( 0 ) = 1, and we 


look for solutions of (4.2) with (/){/3,y) = 02 ( 02 /) and this time 0 < 0. Let us recall 
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the construction of 02 , because it provides some informations needed in the sequel. 
We look for 02 : 1^- —t M such that 02 (^) = 02 (— |^|), 02 ( 0 ) = 0 and satisfying the 
Bessel equation 


0 "(r) + 


N -1 


02 (r) — 02 (r) = 0 for r < 0 . 


(4.6) 


The desired solution of (4.6) is 


i/'zW =(iFi(;t + 1 ; —), 


(4.7) 


where r = N/2 — 1 and oFi{;t + 1; is the generalized hypergeometric function 
dehned as 

on(;r + l;^) y. ( 4 . 8 ) 

^or(r + l + n) n\ 

Indeed 02 ( 0 ) = 1 and, using the fundamental property of the Gamma function 
r(r + 2 ) = (r + 1 ) r(r + 1 ) for every r > — 1 , 

we derive 

2 

V’2(^) = ^^^r00yo^i(;^ + 2;^) (4.9) 


and, thus, 02 ( 0 ) = 0. Moreover, from (4.8) we have that 02 (r) and 02 (r) are dehned 
for all r < 0 , the former being positive and the latter negative. 

By plugging {u,v) = into (4.1) and proceeding as in the pre¬ 

vious case, we obtain the following system: 


-da'^ — dl3‘^ -f ca = /'(O) 

2 I —/id002(0i?) 


—Da^ + ca = 


7 = 


dj3'ilj2{fiR) + z/02(/9i?) 
h 


(4.10) 


ci002(0i?) Z/02(/3i?) 


Observe that, in this case, 7 > 0 without any restriction on 0 < 0. We now extend 
the previously dehned functions a^(c, 0) and q;^(c, 0) to negative values of 0. Solving 
the second equation of (4.10) for a gives 


a 


where 


i(c, /^) = 0 ^ ± -411x2(0)) , 


(4.11) 


d(3'ip2{(3R) + z/02(0i?) ’ 

The function X 2 vanishes at 0, is dehned and negative for every 0 < 0. Moreover 02 
is log-convex (for the proof in the case > 1 see e.g. [12], while the case iV = 1 
follows from direct computation, since in this case 02 (r) = cosh(r)), which implies 
that X 2 is increasing for 0 < 0 . As a consequence, 0 h->■ a(0(c, 0) is decreasing, for 
0 < 0. Finally, it is easy to see that the region of the second quadrant delimited by 
this graph and the 0 -axis invades monotonically such quadrant as c t -f-cx). 
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Now the first equation of (4.10) describes, for c > ckpp, a half-circle (/3, a^{c, /3)), 
(3 <0, which has center (0, and radius 


( \ c|pp 

''W = ' 

in particular it is dehned for 0 > /? > ,d(c) := —p(c). The half-disk bounded by it 
shrinks to its center as c ckpp and invades monotonically the whole of the second 
quadrant as c — >■ -1-cxd. 

In sum, we have dehned the functions for (3 ranging in the whole real line. 

Direct computations show that, for every c > 0, 

limaj(c,/?) = \\mdfja^{c,l3) = 0, (4.12) 


and, for c > ckpp. 


lima^(c,/3) = (d/5) = 0. (4.13) 

As a consequence, the sets S£)(c), Srf(c) dehned by 

^d{c) := {(/?,a^(c,/3)),/3 < /3(c)}, Srf(c) := {(/3, a^(c,/3)),/3 > /3(c)}, (4.14) 


are continuous curves in the (/3, a) upper half-plane, which are moreover diherentiable 
if c 7^ Ckpp- Thanks to the above mentioned monotonicity properties in c of the 
functions and there exists c* such that the curves T^nic) and Tjd{c) do not 
intersect for c < c*, touch for the hrst time, being tangent, for c = c* and are secant 
for c > c* (see Figure]^. The (/3, a) corresponding to the intersection points provide 
the desired plane waves, for any c> c*. 

Actually, we wish to know how many intersection points the curves and 
exactly have for c > c*. The answer is given in the following result, which, apart 
from the importance that it will have in the rest of this work, is interesting because 


it provides information about the number of solutions of systems (4.3) and (4.10) 
using PDF tools, specihcally the comparison principle. We think that this kind of 
technique may be crucial in the study of other reaction-diffusion systems. 


Proposition 4.1. The curves S£)(c) and Srf(c) defined in (4.14) have, fore > c*, at 
most two intersections. Moreover, for c = c*, there is a unique intersection. 


Proof. To prove the hrst part, assume by contradiction that there exist (/3j,aj), i G 
{1,2,3}, with /3i 7 ^ 13j for i ^ j. By the monotonicity in fi of the curves, which 
has been established above, we have that a, 7 ^ aj for i ^ j and we can label the 
intersections so that 

0 < cti < 02 < 03 - (4-15) 


As a consequence, problem (4.1) admits three distinct solutions {ui,Vi), i G {1,2,3} 
of the form (4.2) with {a, (3,^) = {ai, (3i,'yi), where 7 * = 7 (/ 3 j) is given by ( 4.4[) . The 
function {u,v) := {u 2 — ^ 3,^2 — ^^ 3 ) is also a solution to (4.1). Thanks to (4.15), and 
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since 7 i 0 (/ 9 i, y) > 0 in R), there exists L > 0 such that, for every y G -BAf(0, R) 

and t = 0, 

{ui,Vi) > {u,v) for X < —L, ("h, h) < (0,0) for x > L. 

As a consequence, for large fc > 0, we have k{ui,vi){x,y,0) > {u,v){x,y,0) in 
Setting 

k* = mm{k{ui,vi){x,y,0) > {u,v){x,y,0) for {x,y) G f^}, 

k>0 

there exists (xo,|/o) ^ ^ such that 


k*ui{xo, I/O, 0) = u{xo, I/O, 0) 


Now, the conclusion of Proposition 2.1, which obviously holds true for the linearized 


system (4.1), ensures that k*{ui,vi) = {u,v), which is impossible because ui is posi¬ 
tive while u is positive for large negative x and negative for large positive x. 

Passing now to the second part of the statement, assume by contradiction that 
for c = c* the curves intersect in two distinct tangency points. By continuity and 
the monotonicity in /3 and c, each tangency point gives rise to two intersections for 
c> c*, c ^ c*, which is excluded by the hrst part of the Proposition. □ 

According to the position of the hrst intersection of the curves we give the following 
Definition 4.2. Denoting by {(3*, a*) the tangency point between the curves Tjd{c*) 


and Ec;(c*), which is unique by Proposition 4.1 We say that 


c* is of 


type 1 
type 2 
mixed type 


if 13* > 0 ; 
if 13* < 0 ; 
if 13* = 0 ; 


we write c* = 
we write c* = C 2 
we write c* = c* 


(cf. Figure 
(cf. Figure 


B)) 

E)) 


As an important consequence of Proposition 4.1, we can characterize the type of 
c* according to the different parameters of the problem. This will be used in Sections 
and 1 ^ to study some important properties of the asymptotic speed of propagation. 
We start with the case D < 2d and c = ckpp, for which the hrst relation of (4.12) 


and (4.13) gives 


«z)(ckpp,0) = 


Ckpp 


> 


Ckpp 


= a|(cKPP,0). 


(4.16) 


D - 2d 

This entails that c* < ckpp if D < 2d and, therefore, it is of type 1, since ^^(c*) 
is dehned for (3 > (3{c*) > 0. The same conclusions hold true in the case D = 2d, 
where the curve ^^^(ckpp) and the degenerate hyperbola Srf(cKPp) are secant, because 
«d(ckpp, 0 ) = a^(cKPP, 0 ) and 5 / 3 aJ(cKPP, 0 ) = 0 , while c?/ 3 a))'(cKPP,/d) is a negative 
constant for /3 > 0 . 

For D > 2d, instead, if we set 


D 


Cm — 


^JM{D - d) 


Ckpp, 


(4.17) 


thanks to the hrst relations of (4.12) and (4.13), we have that 


Ed(c) n Srf(c) n {/3 = 0 } 7 ^ 0 if and only if c = cm- 


(4.18) 
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Figure 3: Relative position of S£)(c) and Srf(c) as c increases (from left to right): c* = 
(first row) and c* = C 2 (second row). 


We will now show that the value of the second derivatives with respect to (3 of the 
left and right branches of and evaluated at (cm,0), characterizes the type of 
c*. Namely, if 


—2dpR 

Nucm 


lim 

/ 34.0 d/3 


a^{cM,/3) >lim7— {cm,/3) 


/3to df3‘^ 


-2d 


'J' 


Cm Ckpp 


then Sd(cm) and 'L(i(cm) intersect, apart from /3 = 0, for some /3 > 0. Proposition 


4.1 now entails that the curves do not intersect for {(3 < 0} and the monotonicity in 


c of the curves implies that c* is of type 1. Thanks to (4.17), the above relation reads 
as 

D(Nv — pR) > —2dpR. 

Similarly, if 


2dpR d^ 2d 

= te P) > l.m (c„, P) = 


CrPP 


which can be equivalently written as 


D(pR — Nip) > 2dpR, 

then T,d(cm) and T,ci{cm) intersect for /? = 0 and for some (3 < 0 and c* is of type 2. 
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From this discussion, we have 


c = < 


2d ^ Nv 
fxR 


2d ^ Nu 


(4.19) 


and that, when (]* = 0, necessarily c* = c*^ = cm and 

q2 q2 ^ q2 _ 

lim ,3) = lim (c„. /)) = - lim ^aS(cM. /)) = - lim (cm. P). 

(4.20) 

In the case D < 2d, we have seen in (4.16) that c* = c^, and therefore (4.19) is valid 
for every choice of the parameters. 

With the plane wave solutions constructed in this section we are able to give the 


Proof of Theorem 1.1'^. Let c* be the quantity constructed above and {u,v 


corresponding positive solution of the form (4.2) to the linearized problem (4.1). 


By the Fisher-KPP hypothesis (1.2), for /c > 0, k{u,v) is a supersolution to (1.1). 


Moreover, since the initial datum has bounded support, k can be chosen so tha 
{uo{x,y),Voix,y)) < k{u{x,y,0),v{x,y,0)) 


the 


and, since (0,0) is a subsolution of (1.1), we have from the comparison principle of 


Proposition 2.1, that 


(0, 0) < {u{x, y, t),v{x, y, t)) < k{u{x, y, t),v{x, y, t)), for all t > 0. 


Pick c > c*. Since the a in (4.2) is positive, we have that 


0 < sup u{x,y,t) < sup 

—x>ct —x>ct 


as t ^ +CX). 


'\=R 


The limit for v follows analogously. By the symmetry of the system, the same limits 
hold true in the case x > ct. □ 


4.2 Generalized subsolutions and proof of (1.6) 


We construct now some generalized, in the sense of Proposition 2.4, subsolutions 


to (1.1) which are compactly supported and stationary in the moving frame with 


velocity —c along the x-direction, that is, for system (2.1). To do so, we will make 
use of the construction of c* carried out in Section 14.11 and of a variant of Rouche’s 


theorem from complex analysis, whose proof is given in Appendix 

Proposition 4.3. Let c* be the quantity constructed in Section M Then, for every 
c < c*, there exist c G (c, c*) and two functions u G C{dTl), y G C(f2) which are 
nonnegative, compactly supported, invariant by rotations in y and such thate{u,y) ^ 


(0,0) is a generalized stationary subsolution of (2.1) for all e G (0,1]. 
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Proof. Consider the linearization aronnd n = 0 of the stationary version of (2.1), i.e. 


-dAv{x, y) + cv:,{x, y) = /'(O)n(a;, y) for (x, y) eVt 

-DAq^u{x, y) + cu^{x, y) = uv{x, y)-yu{x, y) for (x, y) G dPt 
'[ddnV^x^y) = yu{x,y) — z/n(x,|/) for {x,y) G dfl. 


(4.21) 


Snppose hrst that c* is of either type 1 or 2, that is, c* = c* with j G {1, 2} (cf. Def¬ 


inition 


4.2). The stndy performed in Section 4.1 concerns the existence of solntions 


to (4.21) of the form 


{u,v) = e“*(l,70j(/3|/)). 


(4.22) 


with q ;,7 > 0, /? G M and (pj dehned there. Observe preliminarily that the fnnction 


fjj appearing in the dehnition of in (4.5) or (4.11) is analytic and that the radins 


of convergence of its Maclanrin series is oo. This is a consequence of (4.8) for both 


-02 and ' 01 , because it can be easily seen, following analogous arguments as in Section 
4.1, that "01 (r) = o-^i(; n -|- 1 ; —r^/4). As a consequence, we can use the same power 


series to dehne (j)j{(3y) ■= ^jJj{f^\y\) for complex 0 too, and the search for plane waves 


in the form (4.2), but now with a, 0,7 G C, works exactly as in Section 4.1, and 


leads to the systems (4.3) or (4.10) even in this case. Solving such systems amounts 


to hnding intersections between the curves T^nic) and ^^(c) from (4.14), with now 
and dehned for 0 G C. By the dehnition of c* given in Section 4.1, E£)(c) 


and Sd(c) have real intersections if c > c*, a real tangency point ( 0 *,a*) for c = c* 
satisfying 0* 7 ^ 0, and no real intersections if c < c*. Complex intersections are 
sought for c < c*, c ^ c*. If the tangency is not vertical, then the function 

h{^, t) := a'l{c* + ^, 0* + r) - a^(c* + ^, 0* -1- r), 

(with k, I G —} that have to be chosen accordingly to which branches are tangent 
for c = c*) when restricted to is analytic in a neighborhood of ( 0 , 0 ) and satishes 
the hypothesis of Theorem A.l Indeed, the strict monotonicity in c of T^oic) and 


T^d{c) yields that ||( 0 , 0 ) has a sign; moreover, since h( 0 , ■) does not change sign, the 
hrst nonzero derivative of h with respect to r at ( 0 , 0 ) has even order and opposite 
sign to ||(0, 0). If the tangent at (0*, a*) is vertical, which happens for D = d,'we can 
parameterize, thanks to the implicit function theorem, the curves T^d{c) and Srf(c) 
using a instead of 0 , obtaining two analytic functions whose diherence h satishes the 
hypothesis of Theorem A.l[ 

Consider now the case c* = d, 


i.e. 


when 0* = 0. Observe that, thanks to (4.18), 


we have c* = cm- In this case, and 'Ed are constructed by attaching two diherent 
analytic curves at the tangency point and therefore the above dehned function h is 
no longer analytic around (0, 0). We overcome this difficulty by considering just the 
curves a( 0 (c, 0 ) and af{c,/3) for c ~ c* and, say, 0 > 0 , extending them analytically 
for 0 < 0 (which amounts to taking their even prolongation) and then dehning h as 
their diherence. Observe that, in this case. Theorem |A.l can be applied with n > 2, 
since (4.20) implies that J^h(0,0) = 0. 

Summing up, in all the cases, we obtain a solution of h(^, r) = 0 for any ^ < 0, 
^ ~ 0, and r = r(^) G C. Namely, for c < c* close enough to c*, we have a solution 
0 , 0,7 G C of (4.3) or (4.10), and moreover, by (A.4) (and also (A.5) in the case 
0* = 0), 0 (respectively a when D = d) has nonzero real and imaginary parts. Using 
the hrst equation of the systems, we then get that Imo 7 ^ 0 also ii D ^ d. 
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In the end, we have a complex solntion 
linearity we conclude that the pair 


[u^v) of (4.21) in the form (4.22). By 


{ui,vi) := Re [(M,n)] 


is a real solution of (4.21). 

Of course, the above construction can still be achieved if /'(O) is replaced by 
/'(O) — 6, with 0 < 5 < /'(O). This penalization only affects the quantity c*, which 
is easily seen to converge to the original one as 5 0. Hence, for any c < c*, we 

can construct {ui,Vi) solution to (4.21) for some c G (c, c*) and /'(O) penalized by 
(5 > 0, 5 ~ 0. Let us denote := Rea, a* := Ima and ip{y) := 'j(()j{(3y) = ■y'ijjj{(3\y\). 
Direct computation shows that 


ui{x,y) = e“""'cos(aix), vi{x,y) = e'^’''"\^p{y)\ cos{aiX + ATg^p{y)) 


where Arg denotes the argument of complex numbers. 

It is then apparent that the regions {ui > 0} and {ui > 0} fl dVt are just one 
the translation of the other. Moreover, hxing a connected component E of {ui > 0}, 
there is at most one connected component F of {ui > 0} intersecting E] if there is 
none, we take as F any component of {ui > 0}. Then, calling E := E x [0, +oo) and 
F := F X [0, +cxd), for all £ > 0, the pair (u, y) dehned by 


{u{x,y,t),y{x,y,t)):= 


e{ui{x, y),vi{x, y)) in E x F 
(0,0) otherwise. 


satishes condition (2.2). Moreover, if e is small enough, f{evi) > (/'(O) — and 

therefore {u, y) fulhlls all the requirements of the proposition. □ 


We now have all the elements to give the proof of Theorem 1.1 ii), but, before, we 
recall a result from |1] that will be needed. Actually, the statement in [Ij is related 
to the case of the half-space (i.e. problem (1.4)) but the proof directly adapts to our 
case. Anyway, for the sake of completeness and in order to remedy some typos in |1], 
we will provide it. 

Lemma 4.4 (|1], Lemma 4.1). Let ci ^ C 2 be such that any nonnegative, bounded 
solution {u,v) ^ (0,0) of (1.1) satisfies, fori G {1,2}, 


lim u{x + Cit,y,t) = -, 

t^ + CXD H 

lim v{x + Cit,y,t) = 1, 

i^+cxo 


locally uniformly in {x, y) G dVL, 
locally uniformly in {x,y) G D. 


Then, 


lim sup 

\y\=R 


u{x,y,t) — 

/i 


= 0, lim sup \v{x,y,t) 

4 - 5.+00 cit<x<C2t 

\y\<R 


11 = 0 . 


Proof. Let {u,v) be as in the statement, £x e: G (0, i///i) and consider the solutions 
(M^,n^) of (|1.1 ) starting from {uI,Vq) := (sup m, supn) and starting from 
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(mqjVq), with 0 < Mq < z///i — e, Uq ^ 0, supptiQ C [—1,1] x {||/| = R} and = 0. 
By hypothesis, there exists T > 0 such that, for i G {1, 2}, 



z/ 

e 


v'{cit,y,t) - 1 

sup 

u {cit,y,t) - 

^ 2’ 

sup 

\y\=R 

t>T 

y 

\y\<R 

t>T 



(4.23) 


For the same reason, calling k := max(l, |c 2 — ci|T), there exists T' > 0 such that 


sup 

\x\<k, \y\=R 
t>T' 


u(x + C 2 t,y,t) - 


< £, sup \v{x + C 2 t,y,t) — 1\ < e. (4.24) 

|a:|<A:, \y\<R 
t>T' 

Fix A G [1/2,1] and consider c = (1 — A)ci + Ac 2 and t > 2T'. 


If (1 — X)t < T, then applying (4.24) with x = {c — C 2 )t G [—A;, k] yields 


sup 

\y\=R 


u{ct,y,t) - 


< 


sup \v{ct,y,t) — 1| < £. 
\y\<R 


(4.25) 


If instead (1 — X)t > T, we have, by construction and (4.24) again, 

ul{x,y) < l[_i,i](a;) < u{x + C2Xt,y, Xt) < ul{x,y), for all {x,y) G dkl, 

y) < + C2Xt, Xt) < Vq{x, y), for all {x, y) G 12. 


By the comparison principle, considering the evolution by (1.1) of the above data 
after time (1 — X)t and at a; = Ci(l — A)t, we derive 

M^(ci(l - X)t,y, (1 - X)t) < u{ct,yR) < v}{ci{l - X)t,y, (1 - A)t), 

n^(ci(l - A)f, y, (1 - X)t) < v{ct, y, t) < v^(ci(l - X)t, y, (1 - X)t). 

Property ( |4.23[ ) then allows us to infer (4.25) also in this case. Since A ranges in 
[1/2,1], this shows that, as t — )■ +cxd, {u,v) converges to (u/y,!) uniformly in the 
section of the cylinder with x between t and C 2 t. By exchanging the roles of Ci 
and C 2 , we obtain the uniform convergence in the whole section between cR and C 2 t, 
as desired. □ 


Proof of (1.6). We are going to prove that, for all c G (0,c*), there exists c G 
(c, c*) such that any solution {u,v) to (1.1) with bounded nonnegative initial da¬ 
tum (uo,vo) ^ (0,0) satishes 


lim u(x ± ct, y, t) = —, 

t —>^ + 00 jj^ 

lim v(x ± ct, y, t) = 1, 

t^+oo ^ ^ ’ 


locally uniformly in [x, y) G dfl, 
locally uniformly in {x,y) G 12. 


(4.26) 


Then, (1.6) with c arbitrarily close to c* will follow from Lemma 4.4 applied with 
Cl = —c = —C 2 . Of course, by the symmetry of the problem, it is sufficient to prove 
(4.26) in the case of propagation towards left, i.e., x — ct. To do this, we consider the 
pair 

(u(x, y, t),v{x, y, t)) := {u{x - ct, y, t),v{x - ct, y, t)). 
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This is a solution to (2.1) with initial datum From the comparison principle 

of Proposition 2.1 we have that at, say, t = 1, {u,v) > (0,0). By Proposition 4.3 
there exists c G (c, c*) and a generalized stationary subsolution {u,v), in the sense 
of Proposition 2.4, to ( |2.1[ ) which is rotationally invariant in y and lies below {u,v) 
at f = 1. By Proposition |2.4[ this order is maintained for all later times, and from 
Propositions |3.1 and 3.3| we obtain that 


u 


,1 <\immi{u{x,y,t),v{x,y,t)) 

y J t^+oo 


locally uniformly in hi. The proof of (4.26) in the case x — ct is thereby achieved 
thanks to Corollary 3.4 and Remark 3.5 ii). □ 


5 Limits for small and large diffusions 

In this section we will study how the speed of propagation behaves as a function 
of D, the diffusion on the boundary of the cylinder. For this reason we will denote 
it by c*{D). The next proposition yields the first relation in (1.7), as well as the 


characterization of the limit Cq as the speed of propagation for the semi-degenerate 


problem (2.3). Notice that the latter is formally derived from the system (1.1) by 


letting H J, 0. 

Proposition 5.1. The function D i—>■ c*{D) is increasing and, as D fO, tends to the 


(positive) asymptotic speed of spreading of the semi-degenerate problem (2.3) 


Proof. It is easily seen from (4.5) and (4.11) that the function D i—)■ af){c,/3) is 
increasing if /3 > 0, while D i—)■ af) is decreasing for every /3 G M. Hence, since the 
curves and T,d{c*{D)) are tangent, they will not touch for c = c*{D) and 

D' > D. As a consequence c*{D') > c*{D), which gives the desired monotonicity. 

Let D < d and recall from (4.16) that c*{D) = c\{D) < ckpp in such case. We 
further have that 


c c 

minOnfcTd) = -f- > ^ = maxoTfcI, d). 

It follows that the tangent point between T,d{cI) and Sd(c^) is actually between the 
graphs of •) and a^(c^, •)• As H 0, af) tends locally uniformly in /? G [0,/d) 

to the function 

- ( fp ■= = -yd(3'if[{f3R) 

■ c c{d^ij[lftR) + ’ 

which is increasing in fd and satishes 

«b,o(c, 0) = 0, lim ab,o(c, /^) = +oo- 
0t0 

Moreover, the curve afy^^c,-) increases to -foo as c J, 0 and therefore it does not 
intersect aj"(c, •) for c ~ 0. On the other hand, the curves are secant for c = ckpp 
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and, as a consequence, they are tangent for some value Cq G (0,ckpp). This is the 
desired limit of c*{D), owing to the locally uniform convergence of to a^ Q. 

To see that cq coincides with the asymptotic speed of spreading for the semi¬ 


degenerate problem (2.3), it is sufficient to consider the plane waves 


{% v) = 




and repeat the arguments of Section to construct supersolutions and generalized 
subsolutions for the system. Indeed, plugging (u,v) in (2.3) one is led to hnd re¬ 


spectively real and complex intersections between the curves and fhe limit 
of ab as D I 0 introduced above. This is precisely the way in which cq has been 
dehned. With the super and subsolutions in hand, in order to conclude one has only 


to ensure that the comparison principles hold true for (2.3) and the analogous version 


with an additional transport term in the a:-direction. The weak comparison principle 
can be derived using the standard technique of reducing to strict sub- and superso¬ 
lutions and then repeating the arguments of [31 E], while the strong one is given by 
Proposition 2.5 and Remark 2.6 □ 


We complete this section by studying the behavior of c*{D) as D goes to -|-cx). 


giving the proof of the second relation in (1.7) 


Proposition 5.2. As D increases to -|-oo, c*{D) increases to -|-oo too, and 

c-{D) 


0 < lim 


D—^+oo \/T) 


< OO. 


(6.1) 


the semi-degenerate problem (2.4). 


More precisely, the limit in (5.1) eoineides with the asymptotic speed of spreading of 


Proof. We will prove that both cl{D) and C 2 {D) increase to -fcx) as H —>■ -|-cx) and 
both satisfy ( 5.1[ ). Then, the conclusion of the proposition will follow from (4.19). 

The case of cl{D) follows as in the proof of [THl Proposition 7.2], with little 
modihcations. 

The case of C 2 is more involved, since it requires some asymptotic expansions 
of the generalized hypergeo metric function dehning ' 02 ; for this reason we give the 
details. From the monotonicities in (3 of af, and af we have (see Figure [^E)) 


ad(c2p),0) < lim af^icKD), (3). 

p^—OO 


(5.2) 


To calculate the limit in (5.2), we recall the well known relation between the hyper¬ 


geometric function qFi and the Bessel J function ([lU page 100]) 

1 


Jr{z) = 


oFi{; r -I- 1; — —), 


F(r + 1) V2 

which, owing to (4.7) and ( |4.9 ), entails 

02(r) _ 2(r + 1) oEi(; J^{ir) 


0^(r) 


= i 


Fi(;r + 2;’^) Jr+Sr)' 


24 



























This, together with the following asymptotic expansion for Jr{z) in a neighborhood 
of \z\ = oo (see m § 7.21]) 


Mz) = + " (i)) + 1)^) + O Q) sin [z - (2r + 1)^) , 

leads to 

'ilj 2 {r) , cos (zr - (2r + 1)|) ie 2 r+i( 2 r+i)f 

lim ,, , , = i lim - 7 -^^^ = — lim - „ , = —1. ( 

r-5.-oo 'ljj 2 {r) r-s.-oo gjn yir — (2r + l)f) r ^-00 \ _ g2r+i(2T+l)2 


As a consequence, the function X 2 (/d) in the dehnition (4.11) of tends to —/r as 
[3 —>• — 00 , whence (5.2) reads as 


1 



' 1 - 


'-KPP 

c!/(D) 




'l+4;i 



(5.4) 


If 02 ( 11 ) was bounded as H —)■ + 00 , (5.4) would give a contradiction (recall that 
D I— >■ 02 ( 11 ) is increasing and larger than okpp), so c^iD) tends to +00 and (5.4) can 
be written as 


-KPP 

2d 


+ 0(1) < 


cf{D) 


D 


1 +Wl + 4/i 


D 


cf{D)^ 


(5.5) 


On the other hand, we have (see Figure |^E)) 2 ( 02 ( 71 ), 0) < ck^ 2 ( 02 ( 71 ), 0), which 

gives 


1 1 
Z1 ^ M 


1 - W1 - 


-KPP 


of (71) / ’ 


that is. 


D Ad ^ ’ 


Therefore, taking the limsup as H —)■ +00 in this relation and the liminf in (5.5), we 
obtain 

, ,cf{D) , cf{D) 

0 < liminf ——— < limsup ——— < 00 . 

D^ + 00 U £)^_|_oo D 

It is then natural to perform the change of variables 


0 = 




01 


= cxa/D 


in (4.10), obtaining 


.OK 


+ 0(5 = /(O), 

~2,~~ -fid/3ij'2{/3R) 

—a +001 = 


dj3^/J2{/3R) + i"iIj2{(3R) 
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For /? < 0, the second equation describes the curve S£)(c) introduced in (4.14) with 
D = 1, therefore the a solutions of the system are bounded independently of D. 
Thus, taking the limit as D ^ +cx) the system becomes 


+ ca = 


—a + ca = 


—lJt,dl3'ilj2{l3R) 

d/3‘ip2WR) + ^‘ip2{PR) ’ 


(5.6) 


whose hrst equation describes a parabola. The hrst value of c, denoted by C 2 , for 


in (5.1) 


which the two curves of (5.6) intersect, being tangent, provides us with the limit 


To see that the limit in (5.1) coincides with the asymptotic speed of spreading 


for the semi-degenerate problem (2.4), one repeats the construction of plane wave 


solutions of Section]^ this time for (2.4). This leads exactly to the system (5.6) and 
the analogous one with 'ip 2 replaced by V’l- Finally, the validity of the needed weak 
comparison principle follows as in the discussion at the end of the proof of Proposition 
5.1, while the strong one is again given by Proposition 2.5 and Remark 2.6 □ 


6 Limits for small and large radii 


This section is devoted to the proof of Theorem |1.2| ii), i.e. to the study the behavior 
of c* as a function of R. For this reason, the dependence on R will be pointed out in 
all the quantities introduced in the previous sections. We begin with the limit i? —)■ 0. 


Proof of the first limit in (1.8). Observe preliminarily that, thanks to (4.19), c*{R) = 
cl{R) for i? ~ 0. The function xi = appearing in the dehnition (4.5) of 

a^ converges locally uniformly to 0 as R 0. Hence, the functions a^ converge 
locally uniformly to the constants c/D and 0 respectively. It then follows from the 
geometrical construction of cl, see Figure]^ and in particular from the fact that Srf(c) 
(which does not depend on R) intersects the /3-axis, that cl{R) —)■ 0 as R 0. □ 

We study now the limit of c*{R) as i? —?■ -fcx), and its identihcation with the 
asymptotic speed of spreading in the half-space, i.e., for problem (1.4). This 


speed is derived in [3] in the case = 1, but the arguments easily adapt to higher 
dimensions, providing the same speed (independent of N) along any direction 
parallel to the hyperplane y = 0. For the sake of completeness, we give some details 
on how this can be performed. 


Derivation of the spreading speed for (1.4) 

dV-I 


Let us focus on the hrst component. 


denoting x = {z,x') G M x Plane waves supersolutions can be constructed 

“ignoring” the extra x' variable, i.e., in the form {u,v) = {u{z,t),v{z,y,t)), reducing 
in this way to the situation = 1 considered in [3]. This shows that the asymp¬ 
totic speed of spreading in the ; 2 -direction cannot be larger than the spreading speed 
obtained in |3], that we denote by c^. 

The same argument cannot be used to obtain the lower bound, because the x'- 
independent subsolution would have an unbounded support, and this will not allow 
us to conclude. To obtain a compactly supported subsolution one can proceed as 
follows. 
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In the case D > 2d, consider plane waves of the form 

{u,v) = , 

where /3 < 0, L > 0 and $ is the Dirichlet principal eigenfnnction of —dA in the 
{N — l)-dimensional ball of radius L, i.e., 


—dAd) = A(L)d> in i?7v-i(0,-h), 

$ = 0 on dBj^_i{ 0 , L). 

The quantity 7 ^ is chosen in such a way that v vanishes ai y = L, that is, 7 ^ = 76 ^^^, 
which tends to 0 as L —)■ +cx). It is well known that A(L) —)■ 0 as L —)■ +cx). As 
a consequence, if we plug {u,v) into the linearization of ( 1.4[ ) around (0,0), we are 
reduced to hnding intersections in the (/3, a)-plane between two curves approaching 
locally uniformly, as L —)■ +cxd, respectively 


« d ( c ,/ 3 ) := ^ {^c + ^/- ADxc 


with Xooif^) ■= 


ydfd 


d/3’ 


(6.1) 


and the half-circle M+ 9 /3 h-)■ a^(c,/3), dehned in Section 4.1 above (see page 15). 


These are the same curves obtained in [3], which eventually lead to the dehnition of 
c^. In particular, exactly as in [3], if iA > 2d, we hnd complex solutions a,/3 of the 
system for c < c^, c ~ c^, and L large enough. The set where the real part of the 
associated {u,v) is positive has connected components which are bounded in z, y, 
and also in x' because <h vanishes on L). From this, one eventually gets the 

compactly supported shifting generalized subsolution. 

In the case D < 2d, we need to construct a compactly supported subsolution 
moving with any given speed 0 < c < = ckpp- This is achieved, in a standard way, 

by simply taking a support which does not intersect the boundary with fast diffusion. 
Namely, call Ac(T) the Dirichlet principal eigenvalue of the operator —dA + cd^ in the 
ball Bn+i{0, L), and <I)c the associated positive eigenfunction. This operator can be 
reduced to the self-adjoint one —dA -|- c^/(4d) by multiplying the functions on which 
it acts by This reveals that \c{L) — c^/{Ad) = Ao(T), which tends to 0 as 

L —)■ -|-cx). Hence, since 0 < c < ckpp = 2^Jdf'{^), 


As a consequence, for L large enough, after suitably normalizing and extending it 
by 0 outside i3Ar+i(0, L), we have that 


{u{x,t), v{x, p,t)) = (0, <Se(z + ct,x', y - L - 1)) 


( 6 . 2 ) 


is a generalized subsolution to (1.4). 


The existence of these subsolutions immediately implies that solutions to (1.4) 
cannot spread with a speed slower than c^. It actually implies more: reasoning as 
in Section]^ one can obtain the Liouville-type result for (1.4) and eventually infers 


that is actually the spreading speed for such problem. 

We now go back to our problem, showing the convergence of the propagation 


speed for (1.1) to the one for (1.4) as i? —)■ +oo. 
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Proof of the second limit in (1.8). We start with the case D < 2d. We know from 


(4.19) that c*{R) = cl{R) in this case. We need to show that lim/{_j.+oo cl{R) = ckpp. 


Using the log-concavity of (fi as in Section 4.1, we hnd that the function R i—>• af,{c, ft) 
is decreasing in its domain of dehnition and that R h->■ af){c, ft) is increasing. Since 
do not depend on i?, this implies that is increasing and that the limit for 

R —)■ -l-cxD exists. 

Fix 0 < c < CKpp. For L large enough, the pair {u,v) dehned before by (6.2) is 


a generalized subsolution to (1.4), and its support is contained in a closed cylinder 
with radius L which does not intersect the boundary with fast diffusion. Therefore, 
if R > L, after a suitable translation, we can have this support contained in the 


interior of the cylinder and then get a subsolution to (1.1) as well. As usual, htting 


this subsolution below a given solution to (1.1) at time 1 and using the comparison 


principle, one derives Ci(i?) > c. Since this holds for any 0 < c < ckpp, the proof in 
the case D < 2d is achieved. 

Passing to the case D > 2d, we know from (4.19) that c* = C 2 for large R. More¬ 


over, (5.3) entails that the curve aj introduced in (4.11) converges locally uniformly 
to the curve ag" from (6.1). By continuity, for any given c > c^, the curve a 




IS 


strictly secant to the half-circle af, and hence the same is true for aj, provided R is 


large enough. We deduce from the construction of c* performed in Section 4.1 that 


lim sup c* (R) < c* 

FI —^-|-oo 


On the other hand, if c < c^, c ~ 


, since the curve ag" lies at a positive distance 
from a^, so does afj for large R and, therefore, no intersection can occur, which gives 

lim inf c*(R) > c* 

R^+00 ^ ^ — 00 


and completes the proof. 


□ 


We conclude this section with the proof of the monotonicities and the rest of 


Theorem 1.2 ii). 


Conclusion of the Proof of Theorem i.^ ii). We have seen in the proof of the second 
limit of (1.8) that R h- )■ cl{R) is increasing. Conversely, the log-convexity of V ’2 


implies that R aJ), considered for /3 < 0, is increasing and therefore R h-)■ c^iR) 
decreases. With these informations in hand, to conclude the proof we only need to 


determinate whether c* = cl (type 1) or c* = (type 2) using the conditions (4.19). 
The hrst condition in (4.19) implies that c*{R) = cl{R) if D < 2d, and thus 


c*{R) is increasing in such case. In the case D > 2d, we deduce from (4.19) that 


c*{R) = cl{R) (increasing) if i? < Rm and c*{R) = cl{R) (decreasing) if R > R 


M, 


with Rm given by (1.9). 


It only remains to prove (1.10), which is equivalent to c*{Rm) = cm- Thanks to 


(4.18), it is sufficient to show that ft*{RM) = 0. From (4.19) we have that ft*{R) < 0 
for R > Rm and ft*{R) > 0 for i? < Rm- As a consequence, up to subsequences 
(recall that the region of the {ft, a) plane where the tangency is possible is bounded), 
we have that there exist ft~^ < 0 and ft~ > 0 such that 


lim ft*{R) = ft 


± 


R^Rt 



































We deduce that tangency occur at both /?+ and jS and, therefore, Proposition O 
entails that = f3~ =0, which completes the proof. □ 

Appendix A Existence of high-order complex ze¬ 
ros for a class of holomorphic func¬ 
tions 

In this section we generalize the theorem given in [3l Appendix B], which proves the 
existence of complex roots for functions which are, in a neighborhood of a hxed real 
root, perturbations of polynomials of degree 2 which depend on a parameter. Here 
we consider functions which behave like polynomials of higher (even) degree and our 


result is used in the construction of generalized subsolutions in Section 4^ The proof 
is based on a parametric version of Rouche’s theorem and is similar to the one in [131 
Proof of Proposition 5.1], however we present it here for the sake of completeness and 
since we use some slightly different estimates. The result reads like follows. 


Theorem A.l. Assume h : —)■ 

and that there exists n G N \ {0} such that 

d^h 


h(0,0) = 0, 


dri 


is analytic in a neighborhood of {^,t) = (0,0) 
that 

(0,0) = 0 /orl<j<2u-l, (A.l) 


/ 


( 0 , 0 )) (^( 0 . 0 )) < 0 . 


(A.2) 


Then, the holomorphic extension of h is such that, for ^ < 0, ^ ^ 0, the equation 

h{f,-) = 0 (A.3) 

admits a complex root r = T{f) which satisfies, for some k > 1 independent of f, 

I ReT(,^)| << Imr(,^) </c|^|^. (A.4) 

Moreover, ifn>2, k> 1 can he chosen so that 

< Rer(0. (A.5) 


Proof. Thanks to the analyticity and (A.l) we have that, in a (complex) neighborhood 


of (0, 0), h admits a Taylor expansion of the type 

r) = + hi(f, T)f + 0(r^’^+^) 

where 

Ah 1 A^^h 

hi = 0(jfj + jrj) 


dh, , 

ai = ^(0,0), 


1 d'^'^h, , 

a2n = 77rwT7rv;::(0,0), 


i2n) \ 


(A.6) 


As a consequence, equation (A.3) is equivalent to 

^2(^, r) := + aif = -hi{f, T)f + 0(r2"+^). 
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Thanks to (A.2), h 2 (^, •) = 0 admits, for ^ < 0, 2n complex solutions of the form 


r,' 


= Tj{^) = j e {0,... ,2n — 1}, with K{^) := The distance 

between two roots Tj, ti, j ^ I, satisfies \Tj — ti\ > for some constant -d only- 

depending on n. Take ^ It follows that, for |^| small enough, the unique 

solution of h 2 (^, r) = 0 in the ball of center tq and radius denoted by B{tq, 
is r = tq. Moreover, for r G 9i?(ro, |.^|^) we have 

2n-l 

l^2(^,'r)| = a2n\r-To\ \r-Tj\ > a2n|'^|‘^('dA:(0 “ 1^1'^)^"“^ > (1 + o(l)) 

i=i 

as 0, for some constant /c > 0. On the other hand, again for r G dB{To, 
recalling the last relation of ( A.6[ ), 

r)^ + 0 (t 2"+1)| <0(1^1 + If |s)|{| + 0(|f |‘+i) = 0{|f I'+s) 

as 0. As a consequence, taking (A <) <5 < i, we have that, for |,^| small enough, 
+ 0(r^”+^)| < \h 2 {^,T)\ on dB{To,\^\^). Applying Rouche’s theorem, we 
obtain that ■) has the same number of roots inside B{to, |.^|'^) as h 2 {^, ■), i.e., one. 

We call it Relations (A.4) and (A.5) immediately follow since 

Rer(0 G [ReTo - |^|^Rero + l^l"^)], lmr(0 G [Imro - |^|^Imro + 


and Rero = K{^) cos(^), Imro = A'(0sin(^) 


□ 
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